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Abstract. We consider the Cauchy problem on nonlinear scalar conservation 
laws with a diffusion-type source term related to an index s g R over the 
whole space R" for any spatial dimension n > 1. Hero, the diffusion-typo 
source term behaves as the usual diffusion term over the low frequency domain 
while it admits on the high frequency part a feature of regularity-gain and 
regularity-loss for s < 1 and s > 1, respectively. For all s E M, we not 
only obtain the LP-L'^ time-decay estimates on the linear solution semigroup 
but also establish the global existence and optimal time-decay rates of small- 
amplitude classical solutions to the nonlinear Cauchy problem. In the case of 
regularity-loss, the time- weighted energy method is introduced to overcome the 
weakly dissipative property of the equation. Moreover, the large-time behavior 
of solutions asymptotically tending to the heat diffusion waves is also studied. 
The current results have general applications to several concrete models arising 
from physics. 
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1. Introduction 

1.1. Main results. In this paper, we consider the Cauchy problem on a scalar 
conservation law with a diffusion-type source term, taking the form of 

( dtU + V ■ f{u) ^ APsU, xeR",t>0, 

(1-1) 

Here, u = u{x, t) : K" x (0, oo) ^> M is unknown and uq = uo{x) : M" ^ R is given. 
n > 1 denotes the spatial dimension, /(u) = /2(w), • • • , fn{u)) : M — ^ R" 

is a given smooth flux function. Since only the small amplitude classical solutions 
will be discussed, we suppose without loss of generality that 

(1-2) /j(0) = /j(0)=0, l<j<n. 

Otherwise, one can take change of variables 

i=t, Xj ^ X, + f'j{0)t, 

and denote /(•) by 

fiu) = f{u) - /(O) - f'{0)u, 

so that the form of (|l.ip remains unchanged but (jl.2p still holds for /. In the 
diffusion- type source term APgU, Ps with an index s G M is a pseudo-differential 
operator defined by 

PsU^T-'{[m{0]-'Tu}, 
where through this paper the frequency function m(^) is supposed to be strictly 
positive and continuous in .J G R". Moreover, the following assumptions on m(^) 
which describe its more precise behavior as |^| is close to infinity or zero will also 
be used: 

(A^i) Whenever |^| ^ 1 is close to infinity, m(^) Notice that since m(^) 

is strictly positive and continuous in ^ G M", it is equivalent with the 
assumption that m(^) ~ 1 -I- for ^ G M", that is, there are constants 
TT^i > Tno > such that for any ^ G M", 

(1.3) ™o(l + |e|') <m(0 <mi(l + |e|2). 

(A^2) Whenever |^| ^ 1 is close to zero, m(^) — m(0) vanishes with rate |^|'^ for 
a constant a > 0. Precisely, there is m2 > such that 

|m(C)-m(0)| <m2|Cr 
whenever |^| is small enough. 
As will be seen later on, several mathematical models in different contexts can 
provide some concrete examples related to equation (jl.ip i for different forms of 
m(^) satisfying assumptions (A^i) and (A^2)- 

The aim here is to study the well-posedness and large-time behavior on the 
Cauchy problem p.l|) for arbitrary space dimensions n > 1 and for all s G M in the 
framework of small-amplitude classical solutions. Note that under the assumption 
(A^i), 

(1.4) AP,u ^ l ^^Tu] with 



and also the regularity of solutions is generally determined by their high frequency 
part. Thus, it is convenient to say that equation (jl.ip i is of the regularity-gain type 
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for s < 1 whereas of the regularity- loss type for s > 1. In the critical case when 
s = 1, APsU just behaves like a damping term over the high frequency domain; 
more discussions will be given in the next subsection. In this paper, for all s G M, 
we shall establish under (A^i) the global existence of solutions and the time-decay 
rate of solutions and their derivatives up to some order, where the extra regularity 
on initial data is required in the case of the regularity-loss type. Moreover, we shall 
also prove that the obtained solution for all s € R time-asymptotically tends to 
some heat diffusion wave if (7W2) is further assumed. Precisely, these results are 
given as follows. 

First of all, we are concerned with the global existence and time-decay rate of 
solutions to the Cauchy problem (jl.ip under the assumption (A^i). When s < 1, 
one has 

Theorem 1.1 (case of s < 1). Let n > 1, s < 1, and let N > [ti/2] + 2. Suppose 
the assumption (A^i) holds. There is a constant eg > such that if \\uq\\j^n < eo, 
then the Cauchy problem (jl.lj) admits a unique global solution u{x, t) satisfying 

ueC"^ ([0, 00); H^) n ([0, 00); H^-') 

and 

pt 

\\u{t)\\l«+ ||V(V)^"^ii(r)|| dT<C\\uo\\l^ 
Jo 

for any t > 0. Furthermore, whenever |[ito||_ff"nLi small enough, the obtained 
solution u enjoys the time-decay estimate 

(1.5) ||V^-ii(t)|| <Cho||H«nLi(l+i)"^"' 
for any t>Q, where < k < N . 

Remark 1.1. As far as the time-decay estimate (|1.5|) is concerned, the optimal 
L'^ -decay estimate is obtained for all derivatives up to N order. The smallness 
assumption on ||uo||li is used only for the case of derivatives V'^'u with 1 < fc < N , 
but it is not necessary when k = 0; see Remark \3.1\ for details. For s = I, the 
corresponding result in Theorem \l.l\ has been obtained in [4], where the time-decay 
of the L^-norm for u[x,t) only was considered. 

Although the existence result stated as in Theorem 1 1.11 is more or less standard 
for the case when s < for which equation p.ip 1 gains at least one-order derivative 
regularity, we would use a unified elementary energy method as in [4] to prove the 
whole case of s < 1, which in turn can shed light on some main diflaculties for 
the regularity-loss case of s > 1. On the other hand, it is non trivial to obtain the 
optimal time-decay rate of the highest-order derivative u{t) in L^-norm without 
putting additional regularity on the initial data uq. From the proof later on, this 
is actually based on some interpolation estimate on the V^(V)~*-type derivative 
of u(x,t); see Lemma [3731 

Next, when s > 1, the result similar to that in Theorem 1 1.1 1 is stated as follows. 

Theorem 1.2 (case of s > 1). Let n > 1, s > 1, and let N > Nq, where Nq is 
defined by (j3.24p in terms of n and s. Suppose the assumption {A4i) holds. There 
is a constant ei > such that if \\uo\\h" + ||wo||li ^ ^i; then the Cauchy problem 
(|1.1|) admits a unique global solution u{x,t) satisfying 

(1.6) M e c° ([0, 00); H^) n c^ ([0, 00); H^-^) . 



4 



RENJUN DUAN, LIZHI RUAN, AND CHANGJIANG ZHU 



Furthermore, the obtained solution also enjoys the time-decay estimate 

(1.7) 11^=^(011 <C\\uo\\H-nL^il + t)-'i-^ 

for any t > 0, where < fc < Ni with Ni defined by (j3.25p in terms of n, s and N . 

For this time, let us point out some obvious differences of the above two theorems 
corresponding to the case when s < 1 and s > 1. First, for the global existence of 
solutions, the smallness condition on ||ito||Li is needed in the case when s > 1 but 
not required when s < 1. Second, from the definition p.24p . A'o standing for the 
smoothness degree of initial data uq is strictly larger than [^] + 2, which shows 
that the additional regularity of initial data is necessary for the case of s > 1 so as 
to guarantee the global existence of solutions to the Cauchy problem (11.11) . Last, 
from the definition (|3^ of Ni, 

+ 2<Ni<N. 

This together with (|1.7p imply that under the same condition on initial data, 
the optimal time-decay rate in the case of s > 1 is obtained for u and its lower-order 
derivatives only, but not for all derivatives up to iV-order. Here, the optimal rate 
means that it is equal to the one in the case of the linearized equation; see Lemma 
12.11 We remark that all these differences between Theorem 11.11 and Theorem 11.21 
result from the special feature of the equation for different values of s, that is, 
equation p.ip i is of the regularity-gain type for s < 1 but of the regularity- loss 
type for s > 1. 

Finally, we consider the asymptotic behavior of the obtained-above solutions for 
all s G M. It turns out that for either s < lors> l,a properly-defined heat 
diffusion wave can be regarded as a good time-asymptotic profile. For that, set 
fis = l/m{Oy. Note fig > 0. Define the Green function G = G{x,t) by 

G{x,t) = (4/i,7rt)-*e"^ = 7-i{e-^=l«l'}. 
The result is stated as follows. 

Theorem 1.3 (time-asymptotic behavior). Let n > 2, k >Q be integers and s G M. 
Suppose that both assumptions (A^i) and (A^2) on to(^) hold. Let u = u{x,t) be 
the solution to the Cauchy problem (jl.ip constructed in Theorems \1.1\ or \1.2\ for 
given initial data uq G (M") D L^ (IR"); where N is chosen such that if s < 1 
then N > [^] +2 and if s > 1 then N > Nq. Set u* = u*{x, t) by 

(1.8) u*{x,t) = G{x,t + 1) uo{x)dx. 

Then, whenever ||uo||_ff"nLi is further small enough anduQ G L\ with Jjj„ uq dx = 0, 

(1.9) II V'^ {u - u*) {t)\\ < CWuoWH^nLlPim + fj-'i'^-i -^"{i--> 

for any t > 0, where p{t) = ln(l + t) for n ^ 2 and p{t) = 1 for n > 3, and k is 
chosen in terms of n, s and N such that if s < 1 then < k < N and if s > 1 then 
0<k<N2 with N2 defined in . 

Remark 1.2. From (|1.8|) . u*{x^t) carries the same mass of initial data and it is 
usually called the heat diffusion wave. u*{x,t) is a good time- asymptotic profile of 
solutions to the Cauchy problem (jl.ip in the sense that the rate of the nonlinear 
solution u{x,t) converging to the profile u*{x,t) is strictly larger than one of the 



OPTIMAL DECAY RATES TO SOME CONSERVATION LAWS 



5 



nonlinear solution u{x, t) itself decaying to zero due to a > and hence min{l, a} > 
by the assumption (A^2)- 

1.2. Applications and related literature. In this subsection we review some 
known results related to the general model 

(1.10) dtu + V ■ f{u) ^ AP,u 

proposed in this paper. As will be discussed in the following, this general model 
can be connected to some concrete equation or system by taking different values 
of s, which in fact have been extensively studied even for a lot of different issues 
in individual contents. We remark that the exhaustive literature list is beyond the 
scope of the paper, and thus only some closely related results will be mentioned; 
interested readers can refer to them and references therein. 
When s takes the negative value, similar to (|1.4I) . 

APs ^ A - i-Ay-' 

in the sense of the frequency space. Particularly when s = — 1, it is connected with 
the Cahn-Hilliard (CH) equation describing the phase separation in binary alloys, 
in the form of 

(1.11) dtu + Alp{u) + A'^u ^ 0, 

where (p{-) is a given function. Notice that if if is smooth and <p'(0) < 0, the 
linearized equation of (jl.lip is written as 

dtu = -.^'(O)Am - A^w 

which exactly corresponds to the case of s = — 1 in the linearized level. For the 
related study of the CH equation, we mention [1] and [21] only. 

For the Cauchy problem in general n dimensions, |T| studied the following prob- 
lem: 

' dtu = Au+y ^^^-e^A^u, .TGR",t>0, 

(1.12) I f-^^ d^.d^^u 

. u\t=o = "0, X e M". 

Under the assumptions that fij{u) {i,j ~ 1, 2, • • • , n) are Lipschitz continuous and 
equal to a constant outside a bounded interval in u, uo{x) £ L°°(]R"), e\/uo{x) G 
L°°(R") and under some additional conditions, an L°°(R")-a priori estimate, which 
is independent of e > 0, is obtained for the solution of the Cauchy problem ()1.12p . 
On the other hand, few results on the temporal decay estimates and on the asymp- 
totics of the solution of the corresponding Cauchy problem have been obtained. 
Thus two natural questions are under which conditions on the smooth nonlinear 
functions fij{u) {i,j = 1,2,-- - , n) the corresponding Cauchy problem admits a 
unique global smooth solution u(x,t), and how to get the optimal temporal decay 
estimates and further describe its time-asymptotic profile. 

For the above questions, an answer was given in for a special case of the CH 
equation. Precisely, for (|l.lip . when (p'{0) — 0, [21] proved the global existence of 
the C°° solution u for initial data uq G O L°° with ||uo1|li small enough, and 
further obtained the optimal time-decay estimates 

II^MO|Lp(r„) <C^M(l+i)"*"^^'"^^ t>r>0, fc = 0,l,---. 
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where 1 < p < oo. Moreover, for the asymptotic profile, under some additional 
condition on it was also proved in |21) that 



hm (1 + i)'"' 



0. 



As pointed out in [21] Remark 1.2], the case of </5'(0) < can be considered in 
a similar way, which is related to our results for small-amplitude initial data 
stated before. 

When s = 0, one can simply let be an identity operator and thus (|1.10p is 
equivalent with the usual scalar viscous conservation law 

+ V • ![u) = Au. 

The study of the above equation has a long history and has been much more exten- 
sively investigated from different respects, particularly well-posedness and finite- 
time blow up of solutions to the Cauchy problem even for large initial data or 
general flux function, and also stability of wave patterns such as the smooth shock 
wave, rarefaction wave and contact discontinuity, cf. [5]. 

When s = 1, a typical model related to (|1.10[) takes the form of 

-em' 

1+£2|^|2 



(1.13) dtu + y-f{u)^APlu, APlu := \ ^^^^Tu 



where e > is a parameter. This model was derived in |23| as the corresponding 
extension of the Navier-Stokes equations via the rcgularization of the Chapman- 
Enskog expansion from the Boltzmann equation, which is intended to obtain a 
bounded approximation of the linearized collision operator for both low and high 
frequencies. As pointed out in j23[ 125] . we remark that the right-hand term of 
(|1.13p can be understood in the sense that APlu behaves qualitatively like the 
usual Navier-Stokes viscosity Au at the low frequency, while it essentially acts as 
the damping force — u at the high frequency. The first rigorous mathematical study 
of (|1.13p is given in |25] for the one space dimension n ~ I, where the propagation 
of smoothness of solutions with small initial data, existence of traveling waves, 
existence of entropy solutions with BV initial data and zero relaxation limit as 
e — >■ were considered. 

On the other hand, ()1.13|) with e = 1 can be recovered as in [4] from the following 
hyperbolic-elliptic coupled system 



(1.14) 



dtU + V- ifiu) + q]=0, 
~WV ■q + q + Vu = Q. 



When n = 1 and f{u) = hu^ , p.l4p further reduces to a simplified model system 



(1.15) 



dtu + d^i^u'- + q) ^ 
-dlq + q + dxU = 0. 



The above system, which was first derived in [9], actually arises from the study of 
radiation hydrodynamics. A more general system describing the one-dimensional 
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motion of the radiating fluid takes the form of 
' dtp + dx{pu) = 0, 



(1.16) <^ 



dt(pu) + dx{pu^ +p) = 0, 
dt[pie + ^u^)] + dx[pu{e + ^u^) +pu + q]^0, 

^ -dlq + 3a2g + Aand^O^ = 0, 

where the hydrodynamic functions p > 0, u, p, e and 9 > denotes the mass 
density, velocity, pressure, internal energy and absolute temperature of the fluid, 
respectively, and q is the radiative heat flux, and a > 0, k > are the absorption 
coefficient and the Boltzmann constant, respectively; see [37]. 

For some mathematical results on models (|1.14p and (|1.15p . interested readers 
can refer to [4|. Here, we only mention the study of p.l4p in the case of high space 
dimensions. In fact, the model (|1.14|) over R" for any n > 1 was first proposed in 
[6] by taking the proper approximation of the high dimensional version of (|1.16|1 . 
where the one-dimensional result in [20] was generalized to the case of several 
dimensions and also the existence and uniqueness of global weak, entropy solutions 
to the initial value problem and the two different relaxation limits were studied. 
Recently, the stability and convergence rate of solutions near constant states or 
planar rarefaction waves were obtained in [7] [HI [H] for the case when the spatial 
dimension takes values 2 < ?i < 8 on the basis of L^-energy method. For any 
n > 1, [55] exposed the pointwisc estimate of solutions by using the Green's function 
method, and [3] also developed a refined energy method to consider the stability of 
constants states, smooth planar waves and time-periodic solutions in the presence 
of the time-periodic source. Very recently, by using a time- weighted energy method, 
[22] obtained the global existence and optimal decay estimates of classical solutions 
with small amplitude, and also showed that the solution tends time-asymptotically 
to the linear heat diffusion wave. 

Next, for the general model (|1.10p . we turn to the case when s is strictly larger 
than one. If one takes s = 2 and m(^) = (1 + + 1^1^)^^^- similarly as in [3], the 
corresponding model (jl.lOp can be recovered from the hyperbolic-elliptic system 
with a fourth-order elliptic part 

jdtu + V- [f{u) + 9] = 0, 

\ A^g -W -q + q + Vu^O. 

The above system was considered in |10| for one space dimension n = 1, where 
authors established the global solvability and asymptotic behavior of solutions for 
small initial data in H'^ OL^ with N > 7, and also obtained the optimal time-decay 
rates 

\\d'Mt)\\<cii + t)-i'^ 

for < k < [{N + l)/2] — 2. Later on, [TB] considered the one dimensional version 
of the following more general model 

(dtU + V- [f{u) + q]^0, 

(1.17) 

[{~Ayq + q + Vu^O, 

fors > 2, which is indeed equivalent with ([TTOl with m(f ) = (l-|-|^p)^/^ ~ 1 + 
For the model (|1.17[) when 7i — 1 and s > 2, [18] proved some similar results as in 
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[To] . It should be pointed out that both [10] and [18] also showed that the solution 
approaches the nonlinear heat diffusion wave described by the self-similar solution 
of the viscous Burgers equation as time tends to infinity, and furthermore they 
proposed the method to deal with the time-decay property of solutions for those 
equations of regularity-loss type. 

For the regularity-loss phenomenon, it actually has been observed for some other 
realistic systems of equations. For instance, [5] considered the optimal large-time 
behavior of solutions to the Vlasov-Maxwell-Boltzmann system that describes the 
dynamics of the kinetic plasma, and a similar but more elaborate result was also 
obtained in for the Euler-Maxwell system in the context of fluid plasma. In 
addition, |11) and |12j studied the time-decay property for the dissipative Timo- 
shenko system. All these mentioned systems admit a common structure that they 
are of the regularity-loss type. A typical feature that this kind of regularity-loss 
structure generates is that the linear solution semigroup has an upper-bound es- 
timate by e~P^^^*, where the frequency function p(^) behaves like + ICP)^ 
corresponding to the case of s = 2 that we discuss here. We believe that the current 
developed approach could become a powerful tool to provide much more elaborate 
results than those in [5] [3] . 

Finally, we emphasize that although the general model that we propose here can 
cover several concrete examples in different physical contexts mentioned above, we 
work only in the framework of the small amplitude classical solution, which thus 
makes it possible to found a general theory of the global existence and large-time 
behavior of solutions. On the other hand, it could be very interesting to consider 
some other issues in the current setting of the model, such as the well-posedness of 
the Cauchy problem with large initial data and finite time below up of solutions; 
these are left to our future study. 

1.3. Notations and arrangement of the paper. Through this paper, C de- 
notes a generic positive (generally large) constant and A denotes a generic positive 
(generally small) constant. For an integer m > 0, we use iJ™ denotes the Sobolev 
space i?'"(M") with norm || • Wh'", and set = with norm j| • || when m = 0. 
(•,•) denotes the inner product in . L\ denotes the (1 -I- |x|)-weighted L^(M") 
space with the norm 

= / {i + \x\)\f{x)\dx. 

For simphcity, ||/(-,i)||Lp is denoted by ||/(i)||Lp for 1 < p < oo, especially, by 
!/(<)! I whenp = 2. 

The notion u{^) also denotes the Fourier transform J-u{^) of u{x). For a multi- 
index a = (ai, • • • , ttn), we denote 

The length of a is |q;| = ai + • • • + «„. /3 < a means (3i < ai for all 1 < i < n. In 
addition. (V) is defined in terms of the Fourier transform as 

Finally, for r > 0, we define [?']+ by 




r, if r is an integer, 

[r] + 1, otherwise, 
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where [•] means the integer part of the nonnegative argument. 

The rest of the paper is organized as fonows. In Section [21 we study the time- 
decay property of solutions to the hnearized equation by using the Fourier analysis. 
For all s G M, U'-L'^ type time-decay estimates on the linear solution semigroup are 
obtained. The result implies that over the low frequency part, the linear solution 
semigroup keep the same algebraic time-decay rate for either s < 1 or s > 1, and 
over the high frequency part, it decays with some exponential rate for s < 1 but it 
does so with an algebraic time-rate for s > 1 depending on the regularity degree of 
initial data. 

In Section[3l the global existence and optimal time-decay rates of small-amplitude 
classical solutions to the nonlinear Cauchy problem are established. Depend- 
ing on the parameter s G M, we divide the proof by two cases when s < 1 and 
s > 1. For the case of s > 1 which corresponds to the regularity- loss type, the 
time-weighted energy method as in |101 118] is employed to overcome the weakly 
dissipative property of the equation. The key part of the proof in this case is to 
properly define two integers A^o a-nd A^i such that time-weighted a priori estimates 
could be closed. 

In Section HI based on the U'-L'^ time-decay estimates on the linear solution 
semigroup obtained in Section [21 we consider the large-time behavior of solutions 
asymptotically tending to the linear heat diffusion wave (|1.8p . In the last Section[Sl 
for completeness, we use an appendix to collect several lemmas with proofs, which 
will be frequently used through the paper. 



In this section, we study the time-decay property of solutions to the Cauchy 
problem on the linearized equation 



Here, initial data mq = uq{x) is given. It is easy to see that in terms of the Fourier 
transform in x, the solution to (|2.ip is solved as 



As usual, the solution semigroup e'^^°* associated with the linearized Cauchy prob- 
lem (|2.ip is defined by 



Applying Fourier analysis, we now establish the following U'-L'^ estimate on 
the solution semigroup 6"^^=*, which will play a key role in proving the optimal 
decay estimate on the low-order derivatives of the solutions to the nonlinear Cauchy 
problem pTT]) . 

Lemma 2.1. Let n > 1, k > be integers, 1 < p,r < 2 < q < oo, £ > and s G M. 
Define [•]* by 



2. Decay property of linearized solutions 



(2.1) 




u(e,t) = e '"K)=*uo(e). 
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where [■] means the integer part of the nonnegative argument. Under the assumption 
{A4i) on m{^), the solution semigroup e'^^'* of the Cauchy problem ()2.ip satisfies 
the following time- decay property: 



(i) When s <l, 



(2.3) ||V*^e^^=*7.o||L. <C(l+i)-*^^-'^-^hollL.+Ce-^*||V^-+["(- ^"-u^U^ 

for any t > 0, where C = C{n^ k,P, r, q, s) and X ~ A(s) > are constants indepen- 
dent of Mo; 

(a) When s > 1, 

(2.4) llV^e'^^^'uolU, < C{1 + WuoWlp 

+ C(l + i)-^ ||V'^'+[''+"^^-i)l-Mol|L. 
foranyt>0, where C = C{n,k,p,r,q,£,s) is a constant independent of uq. 
Proof. Take 2 < q < oo and an integer fc > 0. From HausdorfF- Young inequality, 



(2.5) ||V'=e^^=*MolU,(Ms) < C 



< C 



\ite-^'u. 



\ite 



L'!'(|{|<_R) 



Mo 



C 



\ite-^'u. 



hi'{\i\>R) 

:= Is,<+Is.>, 



where ^ + ^ = 1 and i? > is arbitrarily chosen. It can be directly observed that 

f |2 

'/L. satisfies the following estimates: 
(i) For any s £ M, there are constants Ci{s,R) > Cq{s,R) > such that 



co{s,R)\e < 



MO' 



<ci{s,Rm' 



holds true over |^| < R. This follows from the fact that m(^) is strictly positive 
and continuous in ^ £ R". 



(ii) For s < 1, there is C3(s, R) > such that 



inf 



C3is,R). 



\^\>R to(^)'* 

For s > 1, there is a constant 6*4(5, R) = ,-ni{i+i/R.}^^ ^ such that 



MO' 



> Ci{s,Rm 



holds true over |^| > R. These two lower-bound estimates are due to the assumption 
(Xi)ofm(O. 

Now, let us make estimates on Is <: and /s,>- First, for Is.< with s £ R, by 
using the standard way, for instance as in [TB] or [TSl Lemma 3.1], it follows from 
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(i) above that 
(2.6) /,,< = C 



< C 



k-Co(sM)\e,\'^t 



We 



li'{\^\<b:) 



Li'm<R) 

< C{l + t)-^^^--^^-i\\uo\\L^(^R^), 

where I < p < 2. Next, we use (ii) to estimate /«,>■ When s < 1, it follows from 
(ii) above that 



(2.7) = C 



L'''m>R) 



< C'e-'^^("'-")* II Id'' 



Mo 



i«'(l?l>fl) 



1^1 



L-'m>R.) 

where Holder inequality 1/q' = (?■' — q')/r'q' + 1/r' with 1/r' + 1/r = 1 for given 
1 < r < 2, and 5 > is small enough. It is obvious that when r = q = 2 which 
implies r' = = 2, it is straightforward to obtain 

Is.> < Ce-C'3(^.«)* |||C|'=?io|L.'(R„) < Ce-^-^(^'«)*||V^^.o||L-(Es). 

Otherwise, when r 7^ 2 or g 7^ 2 which implies r < q and hence (r' — q')/r'q' = 
1/r - 1/g > 0, it follows from that 



Uo 



Uo 



i'-'(l«l>-R) 



Therefore, we finished estimates on I^.y for s < 1, which together with (j2.6p . give 
(|2.3p after plugging them into (|2.5p . In order to prove (|2.4p . the rest is to estimate 
/<,.> with s > 1. In fact, when s > 1, it follows from (ii) above that 



(2.8) Is^> = C 



\^\^e-—'uo 
< C 



Li'(\i\>R) 



< C sup I -^e 
where £ > is fixed. By using the inequality 

1 ~Ci(s,R) 



(l«l>-R) 
^\Mt^'uo\\r.' 



Li'm>R) 



\i\>R 



sup <^ -jye ' ' ' i«i'<='-'' )■ < C(s, £, R){1 + ty^i^ 



for any t > 0, one further has 

/.,> <c(s,£,i?)(i + i)-^||iei'+'"o|L,.(i^i>^). 

Similarly to control > for s < 1, it also holds that 

B\''''M\L^'m>R) ^ Cin^k,i,R)\\S7'^+^'+"^^-'.^^'uohnu^) 
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for 1 < r < 2. Therefore, one has the estimate on /s,> for s > 1 by 
which once again together with 



, prove (|2.4p after plugging them into (|2.5 



This completes the proof of Lemma 12.1 



□ 



It is not clear whether there is an explicit representation of solutions to the 
Cauchy problem ()2.ip since the explicit form of the pseudo-differential operator 
Ps or the frequency function m(^) is unknown. Thus, it is interesting to find an 
asymptotical profile of solutions which has some relatively simple form. For given 
uq, to say that u'^^^{x,t) is a good asymptotical profile of the solution e^^=*uo, it 
means that the rate of e'^^^^uo converging to u^^^{x,t) is strictly larger than that 
of the solution e^^=*uo itself decaying to zero in the same space. In our case, an 
expected choice is the diffusive wave corresponding to the solution to some heat 
equation. In fact, set 



m{0) = 



Then, u''''y(x,i) = e'^^^'wo satisfies 



(2.9) 

For the time-asymptotical rate of c^^^^uq tending to c'^^^^uq, wc have the following 

Lemma 2.2. Let n > 1, k > be integers, 1 < p,r < 2 < q < oo, £ > and 
s G M. Define [£ + — i)]* as in (|2.2p . Under the assumptions (A^i) and (A^2) 
on m(^), the solution semigroup e'^^=* of the Cauchy problem p.ip asymptotically 
tends to the heat semigroup e^°'^* of the Cauchy problem (\2.9\i with some time rates 
stated as follows. Here, in (i) and (ii), a > is given in the assumption (A^2)- 

(i) When s <1, 



(2.10) llV''(e 



)wo1Il. <c(i + t)-^'^-i) 



" \\uo\\lp 



for any t >0, where C — C{n, k,p, r, q, s) and A — A(s) > are constants indepen- 
dent of Uq ; 

(ii) When s > 1, 

( i — i\—h+£. II II 



(2.11) |[V''(e 



)uo\\l^ <Cil+t)-^^^—^> 



\Uo\\LP 



7'=+[»(7-i)l*, 



+ c{i + ty^i^\\v''+^^+"'^r-^)]'uo\\L- 

foranyt>0, where C = C{n,k,p,r,q,£,s) is a constant independent ofuQ. 
Proof. Take an integer k > 0. Consider 



(2.12) 



AP^t 



,At 



In what follows we estimate the above function over the low frequency part |^| < R 
for some R > 0. One can claim that when i? > is small enough, there are 
constants C = C(s, i?) > 0, A = A(s, R) > such that 



(2.13) 
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holds true over |f | < R. In fact, choose a properly small R > such that the 
assumption (A^2) holds, that is 

|m(0 -m(0)| 
|C[<i? \t.\ 



for constants ct > and TO2 > 0. Notice 



so that for any |^| < R, 



1 



Top 



1 



< 



to(^)'* 771(0)" 

^^^^|777(e) - 777(0) |e-^-«l«l'* 



min [r)T,(f)l 



s+l 



<C(,s,i?)|ere-^l«l *, 

where the strict positivity and continuity of ?77(^) have been used, and Xs^b. is 
defined by 

1 



A. 



R = mm 



Therefore the claim mentioned before follows. 

Now, take 2 < q < oo with 1/q + 1/q' = 1. Similarly as in ()2.5|) . by Hausdorff- 
Young inequality, it follows from (|2.12p and (|2.13p that 



klAPA 



e^='^*)uo||L,(K5) <C 



C 



i''(l«l>-R) 



c 



i'''(ki<-R) 



i''(ICI>-R) 



Here, three terms on the r.h.s. can be estimated in the same way as in the proof 
of Lemma 12.11 Precisely, the estimate on the first term is similar to (|2.6p with k 
replaced by k + a, the second term similar to (|2.7p and the third term similar to 
the combination of p.7p and (|2.8p . Collecting these estimates imply the desired 
inequalities p.lOp and p. lip for s < 1 and s > 1, respectively. For simplicity, all 
the details are omitted. This completes the proof of Lemma [221 Q 



3. Nonlinear Cauchy problem 

In this section, we are concerned with the global existence and time-decay rates 
of small-amplitude classical solutions to the Cauchy problem (|l.ip of the nonlinear 
scalar conservation laws for all 77 > 1 and all s S M under the assumption (A4i) on 
777 (^). Throughout this section, wc always suppose the assumption (A^i) and shall 
not mention it for simplicity. 
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3.1. Uniform-in-time a priori estimates. In this subsection, by using the en- 
ergy method as in [4], we obtain some uniform-in-time a priori estimates on the 
solution u{x, t) to tlic Cauchy problem for all n > 1 and s G M. Recall the 

equation 

(3.1) dtu + V ■ f{u) ^ APsU, xeW\t>0. 

In what follows, u[x,t) is supposed to be smooth in x, t and satisfy the above 
equation over < i < T for some < T < cx). 

There are two steps in uniform-in-time a priori estimates. The first step is to 
estimate the zero-order term, and the second step is to consider the energy estimate 
of the derivatives on the basis of Lemma 15.11 when treating the nonlinear term. 

Lemma 3.1. There is \ > such that 

(3.2) |h(t)f + A||V(V)-M^)f <0 
for any < t < T . 

Proof. As in [U Lemma 2.1], the zero-order energy estimate on p.l[) gives 

Id " 

(3.3) 2 ""^^^ + E^/^H-. ' = 

j=i 



By Plancherel theorem. 



{AP,u,u) = {AP,u,u) = - I ^il-lupde 



2 



Here, using (|1.3p . 



\u\ d^<-mm^—,— } . ^ I.I2N.J "! 



Therefore, 



{APsU,u) < -mm{ \\V {V)-'u{t)f 



'0 '"1 

On the other hand, from integration by part. 



{.f]{u)x,,u) = - I II fj{v)dT]> dx = 

for each I < j < n. Plugging the above two estimates into p.3p gives the desired 
inequality (|3T2|) . □ 



Lemma 3.2. Let k >1. There are A > 0, C such that 

d 
di 

for anyO<t<T, where = ^ = ^iia^T.-a^! f""^ " = (ai,a2,- • • ,«»)• 



(3-4) E c^liaMOf + A||vi+'^(v)-Mt)f <C||Vu(t)|U.,|lvMi)f 



|ct| — 

^fc _ fc! _ A;! 
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Proof. Take fc > 1. Similarly as in the /c-order energy estimate on p.ip gives 
1 d " 

\a.\—k J — 1 |a| — A: |a|— 

Similarly before, it follows from Plancherel theorem that 



|a|— IckI— ^ 



\(:\2+2k 



where we used the identity (e? + C| + • • • + f^)'' = E C'^Ci^"'C2"' • • • ■ Then, 

|a|=fe 

from (jl.3p . one has 

(3-5) ^1 E ^'ll^MOf + min|i^,^|||Vi+^V)-^.(t)f 

n 

= E E '-'alj-.a, 
3 = 1 \a\=k 

where for 1 < j < n and \a\ = A:, /j_q, is given by 

(3.6) 4 „ = {d,^d^f,{u), -d'^u) - 

+ ( E CidPr^{u)d'^-Pu,^,-d^vS := + 
\|^l=i / 

It follows from integration by part that 

(3.7) Ilc.-\l {d''ufdx<C\\Vu{t)\\^^ P''u{t)f. 



Next, we estimate ^ for a fixed j with 1 < j < n and |a| = /c > 1. In fact, since 
|/3| > 1, the term f'j{u)d°'~^Uxj can be written as 

with Ci = (0, • • ■ , 0, 1, 0, • • • , 0), 1/3 — Cil + |a — /?| = fc — 1. Therefore, applying 
Lemma |5. II with p = r = 2,q = oo and Corollary 15. 11 one has 

pP-^' {f','{u)u,)d'--"u,^{t)\\ <C\\f'^'{u)u,^{t)\\^^ ||V^-iu,^.(t)|| 

+c||«,,(t)||^„ IIV^-i (/i'("K.)(*)|| 

(3.8) < C||Vu(<)||i^ ||v''u(t)|| . 
Thus, by Holder inequality and (|3.8p . 

(3.9) < C||9"w(t)|j ||9'^/j(w)9"-^zi,,.(0|| < C||V«(t)||^^ ||V^^.(t)f . 
Substituting (P^ . (|X7)) and ((^ into (|53)) gives the desired estimate □ 

3.2. Existence and time-decay rate for s < 1. In this subsection we consider 
the global existence and time-decay rates of solutions to the Cauchy problem (jl.ip 
for the case when s < 1. The main goal is to prove Theorem ll.il 
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3.2.1. Global existence. As the first step, we devote ourselves to the 

Proof of global existence in Theorem ll.lt This follows from the local existence 
and uniform-in-timc a priori estimates as well as the continuity argument. The proof 
of the local existence is standard, for instance, cf. [T7l[T9|, and is thus omitted for 
simplicity. It suffices to consider the uniform-in-time a priori estimates for a smooth 
solution u{x,t) to the Cauchy problem p.l|) over < t < T for some < T < oo. 
From Lemma [3Tl and Lemma [321 by taking the summation of p.2|) and p.4p with 
1 < fc < A^, one has 

, N N 

E C^alia-^^Wf + AEll^'^'^)"Mi)f <C||V«(t)|U.||Vu(t)I|^„_.. 

k=0 \a\=k k=0 

Since 

N N 

^IVI^+^^XV)-" = |V|(V)-^^|V|'^ - |V|(V)-"(V)^ = |V|(V)^-^ 

k=0 k=0 

it further follows that 

(3-10) |E E C^al|5"^.Wf + A||V(V)^-^^(t)f <C||Vu(i)llL^HV.^(t)||?,«-. 

A;=0 \a\=k 

for any < t < T. Let us suppose the following smallness a priori assumption 

(3.11) sup \\u{t)\\HN < 6 

0<t<T 

for a constant S > small enough. Due to > [n/2] + 2, 

sup ||Vw(t)||L=o < C sup \\u{t)\\HN <C5. 
o<t<T a<t<T 

Then, for the case when s < 1, the right-hand term of (|3.10l) is bounded by 

C<5||V(V)^-iu(t)f < C6\\\7{Vf-^u{t)f, 

which from (j3.10p . implies 

, N 

rf^E E C^PMOf + A||V(V)^-Mi)f <0 

fc=0 |a| = fc 

since (5 > is small enough. By noticing 

N 

E c'jd-u{t)r^\\u{t)\\%., 

k=0 \a\ = k 

one can take further time integration over [0, t] to obtain 

(3.12) \Ht)\\l. + f ||V(V)^-^ii(r)f dr < C\\u4l. 

Jo 

for any < t < T. Therefore, in the standard way, as long as ||uo||ff" is sufficiently 
small, the above uniform-in-time a priori estimate obtained under the assumption 
p. lip implies the global existence of solutions by combining the local existence and 
uniqueness. This completes the proof of global existence in Theorem ll.il □ 
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3.2.2. Optimal time-decay rates. In order to finish the proof of Theorem 11.11 the 
rest is to prove (|1.5|) for time-decay rates of the obtained solution u{x, t). For that, 
define 

N 

f°P(i)= sup ^(l + r)t+'=||vMT)f. 

One can prove that £'^{t) is bounded uniformly in time if Huolli/^n-Li is small 
enough. In fact, let us begin with 

(3-13) I J2 C^l|5%Wf + A||Vi+^-(V)-^u(t)f <Cx.>i||Vu(i)||L~l|V^«(t)f 

\a\—k 



for any t > 0, where < k < N. This follows from Lemma [3.11 and Lemma [3.21 
Notice that the right-hand term of p.l3p vanishes when k ~ 0. Let a > be a 
constant to be chosen later. The time-weighted integration of p.l3p gives 

(3.14) il + t)''+^\\W^u{t)f + f {l+T)''+''\\V'+''{Vyu{T)fdT 

Jo 

<C||V'=Mof + C(a + fc) / (l+Tr+^-i||vMT)frfT 

Jo 

+ Cxk>i f (l + Tr+'=||V«(r)!Uo.|lV'=«(r)|l2dr. 
Jo 

To control the second term on the r.h.s. of p.l4p . consider 

llvMOf - iivi+'=(v)-i^t)|p + \\vHy)-'n{t)r- 

< ||V^+i(V)-^u(t)f + ||V'=(V)-^u(i)f 

since s < 1. Corresponding to both parts in the sum above, one need to estimate 
the following two time integrations: 

h=C{a + k) Al + r)°+*'-i||V'^+i(V)-MT)fdT, 
Jo 

h=C{a + k) [ (l + r)"+*'-i||V'=(V)-^?/(T)fdT. 
Jo 

Let a > 1. /i is estimated by 

Ii<S f (l + T)'^+*||V^+i(V)-'u(r)||2fir + C5 / || V^+i(V)-"u(T)f dr 
Jo Jo 

<S f\l + rr+'' II V^+i {V)-'u{t) f rfr + CsWuoWj,. , 



where i5 > is an arbitrarily small constant, the Young inequality with + 
= 1 was used, and the inequality 

E / \\^'+'M-'u{r)rdT<C\\uo\\l. 

k=0''° 



due to (|3.12p was also used. For the estimate on I2, the following interpolation 
inequality similar to [H Lemma 2.4] is needed. 
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Lemma 3.3. For any k > 0, one has 

(3.15) ||V'^(V)^"u(t)|| < C||V'=+i(V)^"u(t)||^™||w(t)||£™. 

Proof. Set 

A=||V^-+i(V)-Mi)ll, B^\\u{t)h^. 

Then, 



5 



1 r |t|2fc+2 /• |c|2fc 



i?2 



< 4^ + Ci?2ii>"+2'=, 

where R> can be arbitrary. By taking i? > such that ^ = CB^R^'^'^'', that is 
it follows that 

which is equivalent with p.lSp . □ 
By the method of [TH [TSl [H] , similar to [H Lemma 2.3], it can be proved that 
(3.16) \Ht)\\L^ < holUi 

for any t > 0. Now, by plugging p.isp into I2 and using p.l6p . one has 



2(Ti + 2fc) 



I2<C / (l + T)'^+'^-i||V'=+i(V)-^u(r)||^T5ET4||ii(r)||£t^'"+^dr 
Jo 

<5 [ (l + T)"+^-||V'=+i(V)-^ii(T)f dr 
Jo 

+Cs f\l + rr+''^\\u{r)\\l^dr 



<s / {i+Tr+'\\v'+'{v)-^uiT)rdT + cs\\uo\\iii+tr--, 

Jo 

where the Young inequality with ^"^^^^ + ^^2k+2 ~ ^ used and a > ^ is 
assumed. Collecting these estimates on Ii and I2 above, the second term on the 
r.h.s. of (|3.14|) is bounded by 

(3.17) Cia + k) f\l + T)''+''-^\\V^u{T)\\^dT 
Jo 

<2S f (l + r)°+'=||V'=+i(V)-^u(r)||2dr 
"'0 

+ Cs\\uo\fH. +Cs\\uo\\liil + tr-9 
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for an arbitrarily small constant S > and a constant a > niax{l,|^}. Moreover, 
the third term on the r.h.s. of p.l4p is bounded by 

(3.18) Cxk>i f\l + Tr+''\\Vu{T)\\L^\\V''uiT)\\^dT 
Jo 

<CXk>l f{l+TY+'^-^\\V'^u{T)fdT^£f{t) 



<Cxk>i (l + r)«+^-i|lV^w(r)fdryfTO: 
where ^^-ji > 1 was used and due to TV > [n/2] + 2, 



\\Vuir)h^<C{l+r)-'^,/£^) 



for any < t < t. Notice that the time integration term on the r.h.s. of (|3.18p can 
again be estimated as in (j3.17p . Thus, under the assumption 

(3.19) supf°P(t)<l, 

t>Q 

by plugging (|3.17p and (|3.18p into (|3.14p and choosing a properly small constant 
(5 > 0, one has 



{l + ty+^\\V^u{t)\\^ + / (l + T)''+^-||V^+'^(V)-"w(T)fdr 
Jo 

<C\\uo\\l.+C{l + tr-'^\\uo\\h 
for any t > 0, where < k < N. Since a > ^, it follows that 

(3.20) sup£^P(t)<C||uofH«nLi 

t>o 

under the assumption (j3.19p . Therefore, by the continuity argument, as long as 
ll^o||_ff«nLi is sufficiently small, p.20p holds true. Then, (|1.5p for time-decay rates 
of the solution u{x,t) follows by the uniform-in-time boundedness of £'^(t). This 
completes the proof of Theorem 11.11 □ 

Remark 3.1. Notice that the right-hand third term of p.l4p can be neglected when 
k = 0. Then, (j3.14p together with (j3.17p imply 



il + tr\\u{t)r+ / (l+r)«||V(V)-Mr)||^dr 
Jo 

<C\\uo\\^ + C\\uo\\l,il + tr-^, 

which further yields 

\\u{t)\\<C\\uo\\L^-nL^{l+t)-^ 

for any t > 0. Therefore, for the above time- decay estimate, the smallness assump- 
tion of ||tio||Li can be removed. This is also consistent with the result in [4] for the 
case when s — 1. 
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3.3. Existence and time-decay rate for s > 1. In this subsection, we consider 
tiie global existence and time-decay rates of solutions to the Cauchy problem 
for the case of s > 1. The main goal of this subsection is to prove Theorem 11.21 
Similarly before, this follows from the local existence and some uniform-in-time a 
priori estimates with the help of the continuity argument. Once again, the proof of 
the local existence is standard and thus omitted for simplicity. In what follows, we 
only consider some uniform-in-time a priori estimates on u(x,t) which is supposed 
to be smooth in a;, t and satisfy equation p.ip over < t < T for some < T < oo. 

In order to state a priori estimates, let us define three integers TVq, N and A^i in 
turn in terms of n > 1 and s > 1 as follows. Notice that n is an integer and s £ M 
is a real number and hence might not be an integer. For any integer fc > 0, denote 



(3.21) 



(3.22) 



hin,s,k) = k+[{- + k)is-l)] + , 



4(n,s,fc) = fc + 1 + 



[(2±3+fc)(s-l)] + , ifn=l, 
[{n±l + k)(s - ifn = 2, 

[(f + fc)(s-l)]+, ifn>3. 



(3.23) 4(5,^) = 2+[2(.s-l)] + fc[.s] + . 

Then, iVo is defined by 

[[^]>[t]+2, fc>^i(n,s,[f]+2), 

fc>^2(n,s,[f]-K2), and A:>4(s,[t]+2) 
N is arbitrarily chosen such that N > Nq, and finally A^i is defined by 

[|] + 2<A:< [^], £i(n,s,fc) <iV, | 

h{n,s,k)<N, and €3(.s,fc)<iV J 



(3.24) 



(3.25) 



No = MikeZ 



Ni = sup < k e 



Now, let us also define some temporal time- weighted functionals £N{t), 2?Af(t), 
f°P(i) andAf,(i) (1 = 0,1) by 



(3.26) 



(3.27) 

and 
(3.28) 



sup (l + r)*^-^||V'=«(r)f^„_.„^, 
t 

J2 I (l+7-)^-5||V'=(V)^-'=W+u(T)||dT, 
i.-n "'O 



fc=0 



fc=0 



0<T<t 



M,[t) = sup (1 + r) — ||V'u(r)||Lo 

0<T<t 



0,1. 



The key point is to prove that all the above functionals are bounded uniformly in 
time if llMoHfl-Npi^i is small enough. In fact, the uniform-in-time a priori estimates 
on these temporal functionals can be obtained in the following 



OPTIMAL DECAY RATES TO SOME CONSERVATION LAWS 



21 



Lemma 3.4. Let n > 1 and s > 1. Let Nq, N and Ni be defined as before. Then, 
one has 

(3.29) £N{t) + VN{t) < C\\uo\\In + CMi{t)VN{t), 
and 

(3.30) f°P(t) < C\\u4l.^^, + C£°^^{tf + Chh{tf£N{t) 
for any < t <T. 

Proof. Wc first prove p.29p . It is equivalent to prove that for any < t < T, 



(3.31) (l + t)^-5 ||v''^u(i)||^„_.[.j^ + / (l + T)'=-i V'^(V)^-'^'W+j/(r) 



dr 



dr 



< CWuoWIn + CMi{t)VNit) 

for all < k < [j^]. This can be done by induction on fc. In fact, similar to obtain 
(|3.10p from p.ip for the case of s < 1, it also holds true for the case of s > 1 that 



N 



dt 



^ ^ c^liaMOf + A||v(v)^-Mt)f <c||v^.(^)|U»||v^.(^)||^„_.. 



fc=0 \a\=k 



Multiplying the above inequality by (1 +t) 2 and then taking integration in t gives 



{l + t)---\\u{t)\\l.+ / (l + r)-^||u(r)||l,Wr 



+ / (l + r)-3||V(V)^^"^w(r)||^dr 



< C ||wo||^,« +C (1 + t)-^ ||Vu(r)||^^ l|Vu(r)||^„_, dr, 
Jo 

where by using the definitions of Mi(t), I?Ar(i) and the fact that n > 1, the last 
time-integration term is bounded by 



C / (l + r)-5l|Vw(r)l|^„l|Vw(r)|i^„_,dT 



<CKUt) / (l + r)-— |h.(T)||^„dr 
Jo 

<C||lio||H«+CMi(t)2?Ar(i). 

Therefore, p.3ip with fc = follows due to s < [s] + . Next, suppose that p.3ip is 
true for fc — 1 with 1 < fc < [p^]- From p.l3p which actually also holds for the 
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case of s > 1 , after taking summation from fc to fc + — fc [s] -|- , one has 

fc+7V-fc[s]+ fc+Af-fc[s] + 



(3.32) 



dt 



i—k \a\—i 



i=k 



k+N-k[s] + 

<c\\vu{t)h^ J2 IIVMOII' 

i—k 



Notice 



A:+JV-fe[s]+ N-k[s] + 

i—k i—0 

which due to s < + , implies 

k+N-kls] + 

Y l|V^+'(V)-''u(i)f > A||V'=+^(V)^-('=+i)W+M(t)|l2. 

i—k 

Hence, multiplying (|3.32p by (1 + t)''^^ and then taking integration in t yields 
(3.33) 

{l + tf-^ ||V'=u(t)||^«_.w+ + / (l+T)^-^||V'=+i(V)^-(^+i)W+w(r)fdr 



+ C (I + t)*^-! V'^(V)^-'^'W+w(t) 



dT 



+ C (l + r)*^-5 |jv?.(T)||^^ V'^-(V)^-'^W + ii(r) 
Jo 

Here, the third term on the r.h.s. of p.33p is bounded by 



dr. 



CMi{t) I {l+rf-"^ V'^'(V)^^'^W + ii(r) 



dT 



<CMi{t) (l + r)'^'~i V'^(V)^-'''W+u(t) 



dr 

< CMi{t)VN{t), 

again from using the definitions of Mi{t)^VM{t) and the fact that n > 1, while 
for the right-hand second term of p.33p . by induction assumption for fc — 1, it is 
bounded by 

C\\uofH^ +CAh{t)VN{t). 

Therefore, from p.33p as well as the induction assumption for fc — 1, p.3ip is also 
true for k. Then, by induction on fc, p.3ip holds for all < fc < [p^]- This proves 

Next, to prove p.30|) . we rewrite (jl.ip as a mild form by Duhamcl's principle. 



u{x,t) = 6^-^^*1*0 - V / 
.7 = 1 



(3.34) 



„AP,(t-r 



)(ff,(«),J(T)dT 
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Here, gj{u) = ,fj{u) — fj{0) — fj{0)u is set for I < j < n. Notice that for each 
J, gj{u) = 0{u^) by the assumption (|1.2p . In what follows, fix an integer k with 
< A: < A^i. Applying V*^ to p.34p and taking norm, one has 



(3.35) II V^u^ll < C II V'=e^^=*Mo|| +CY. j V'^-+ie^^^(*-^)5,(u)(r) 



dT 



dT := /3 4-/4 + /5- 



it follows that 



For h, by applying ([231) with p = 1, r = = 2 and = f + f . 

(3.36) h < C{1 + ty'i-iiWuoWL^ + || V^i("'^^^)^.o||) < C{1 + t)-i~i\\uo\\H-nL^, 

where the definitions p.2ip and p.25p for £i{n,s,k) and A^i were used. For I4, 
one can apply p.4p with k replaced by fc + 1 and with p = l,r = q = 2 and 
2(^s^-i) ~ ""^I^""* + |i where 7(n) = 3 for 71 = 1, 1 for n = 2 and for n > 3, so that 
it follows 

(3.37) h<C [\l + t- T)-t-^ \\g{u){T)\\^, dr 



+ C / '(l + f~r)-^^-f ||V'=+i+[('^+^)(^-i)l+g(^.)(r)||dr 
Jo 

= /4,1 + /4,2- 

Here and hereafter we used g{u) to denote gj{u) for 1 < j < n without loss of 
generality. For the term /4.1, it is easy to see 

ll5(«)(r)IL. < ClKr)f < C£:^P(i)(l +r)-t, 

which implies 

t_ 

(3.38) h.i < C£°^P^{t) /'(l+t-r)-T-^(l + r)-tdr 

Jo 

r(l + t)-i-t, (n = l)j 

<Cf^P(t) X J (l + t)-i-*ln(l+0, (n = 2) > <C£°^^^it)il + t)-^-^ 
[{l + tr^-^, (n>3), 

for all n > 1. 

On the other hand, for the term 74^2, by applying Lemma |5. 11 one has 

(3.39) ha<C r {l + t-T)-''^ 



uiT)\\L^\\V'-^"''^''\L{T)\\dT 
\l+T)-^\\u{T)\\HNdT 



< CMo{t)^/£N{t){l + ty'i-^, 
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where the definitions ([3:281) . ([3:261) and ([3:221) for Mo(t), £N{t) and £2(71, s,k) were 
used and we also used 



(3.40) / {l + T)-^+TdT <C{l + t) — 

Jo 

for all ri > 1. Then, combining (|5351) . ([X^ with ([5371) gives 
(3.41) 



h < C{1 + tr'i-^ [f^P(t) + A/o(t)v/MO 

For the term I^, one can apply (|2.4p with p=:l,7' = g = 2 and ^^^^^^.^ 
a constant (5 > small enough to be chosen later, so that 



1 + (5 for 



(3.42) l5<C [ il + t-T)-^-i \\V''giu)iT)\\^,dT 

+ C j^l + t - r)-(i+^') II Vi+^-+[2(i+*)(«-i)l + g(^)(T)||dT := 7,5,1 + 15,2. 
Here, for Is i, since 

||V^5(u)(r)||^, <C||u(r)||||vMT)|| 



by Lemma 15.11 one has 



(3.43) /5,i<C£^P(t) l{l + t-T)---^-{\ + T)----dT 

({l + tr\-i, (n = l)j 

<C£:^P(t) X J (l + i)-i-lln(l+0, {n = 2) > < Cf^^ (i)(l + i)"^"* 
i(l + t)-t-t, (n>3). 

for all n > 1. 

To estimate 15,2, notice that from the definitions p.23p and p.25p of ^3(5, k) and 
A'^i, one can take 5 > Q small enough such that 

1 + [2(1 + 5){s - 1)]+ + k[s\+ < £3(3, k)<N 
for any < fc < A^i . Then, it follows from Lemma 15.11 that 

||yi+fe+[2(i+A-)(.-i)]+^(^)(^)|j < c\\u{t)\\l^ ||Vi+^+l2(i+*)(--i)l+u(r)|| 
<CMo(0(l + T)-t||v'^7.(T)||^„-.w+ 

for any < r < t. This further implies 



h,2<CMoit)^8^ / (l+t-r)-(i+*)(l + r)-t--+4dr 



<CMn{t)^£^{l + t)-^-^-^ / (l + T)-(i+*)dr 

<CAUt)^£j,{t){l+t)- 
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where n >1 was used. By plugging estimates on 1^,1 and I5.2 into p.42|) . one has 



(3.44) h<C [5^P(t) + Mo(t)v/Mi)J (1 

Thus, together with ([535]) . ([XiT]) and yields that for any < i < T, 



with < fc < A^i . Since the right-hand term of the above estimate is nondecreasing 
in t, it further holds that 



for any < t < T. This proves p.30p and hence completes the proof of Lemma 

inn □ 

Proof of Theorem ll.2t As mentioned at the beginning of this subsection, it 
suffices to consider the uniform-in-time a priori estimates on the smooth solution 
u{x,t) to the Cauchy problem (|l.ip over < t < T for < T < 00. Due to the 
definition p.25p of A^i, iVi > [^] + 2 holds true. By using the Sobolev inequality 
as in [26l Proposition 3.8], it follows that for any < t < T, 



Mo{t) + Ah{t)<C^/£^). 
Then, from Lemma [SH ([XM]) and ([330]) imply 

£Nit) + VN{t) < Cl|uo||^« + C^£^)VN{t), 

and 

respectively. By setting 

X{t) ^ £N{t) + VN{t) + S°^^{t), 
the above two inequalities further lead to 

X{t) < C (hollff^nLi + X{t)i + X{tf) 

for any < t < T. From the continuity argument, it is easy to see that X{t) is 
bounded uniformly in time under the assumption that ||ito||ffNp^i is small enough. 
Therefore, the global existence of solutions as in (|1.6p follows by the standard way 
and also the optimal time-decay estimate (|1.7p results from the definition (|3.27p of 
£'^^ (t) . The proof of Theorem 11.21 is complete. □ 

Up to now, we have obtained the existence and optimal decay rates of the global 
solutions to the Cauchy problem (|l.ip for all spatial dimensions n > 1 and for all 
s G M. 
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4. Large-time asymptotic behavior 

In this section wc shall prove Theorem 11.31 on the large-time behavior of the 
obtained solutions. For that purpose, we divide the proof by several steps in order 
to time-asymptotically approximate the solution to the Cauchy problem p.ip . 

First of all, we prove that the solution to the nonlinear Cauchy problem (jl.ip 
can be approximated by the one to the corresponding linearized problem at infinite 
time. For given uq = uo{x), let us define u = e^^^^uo to be the solution to the 
linearized Cauchy problem corresponding to (jl.ip by 

dtu - APsU = 0, X e R", t > 0, 
u\t=o = uo, X e M". 

Then, one has the following two lemmas which correspond to the case when s < 1 
and s > 1, respectively. 

Lemma 4.1 (case for s < 1). Let n > 2, s < 1, and N > [2^] + 2. Suppose that 
||mo||//np^i is sufficiently small, and u{x,t) is a solution to the Cauchy problem 
(jl.ip obtained in Theorem ll.li Then, for any t > 0, 

(4.1) \\y''{u~e^''^'uo){t)\\<Cpm+t)-i-'^, 
where < k < N , and pit) = ln(l + t) for n = 2 and pit) ~ 1 for n > 3. 
Proof Take < fc < iV. Similar to obtain ((335|) . it follows from (|334)) that 

(4.2) {uit) ~ e^P^'uo)\\ <li + I'„ 

where I[, I'^ have the same definitions as 1^, in (|3.35p . In what follows we would 
improve the previous estimates on I4, with the help of the obtained time-decay 
estimate (|1.5p . For 14, one can apply (|2.3p with k replaced by A: -|- 1 and with 

p=l,r = g = 2, so that 

(4.3) 4<C f\l + t~T)-^-'^\\g{u){T)\\^,dT 

Jo 

+ C r e-^(*--)||V^+ig(^.)(T)|Idr 1 + 4 2, 
Jo 

where for the term I'^ from Lemma [5T] and (|1.5p . 

/•I 

(4.4) Iii<C (l + t-T)-t-^ ||u(r)f dr 

Jo 

t 

< C /'(l + t-T)-T-^(l + r)^tdT < Cp(t)(l + i)"T-^, 
Jo 

and for the term 2 , by applying Lemma 15.11 and (|1.5p , one has 

(4.5) /^2<C /%-^(*-^)||u(r)||io.||V'=+iw(T)||dT<C(l+t)-T-^. 

Jo 

Then, plugging and into gives 

(4.6) /4 < + 
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For the term /g, one can apply p.3|) with p = 1, r = q = 2, so that 
(4.7) I^<C f (l + t-r)-5-^||v'=.9(«)(r)||^,dT 

+ C f e-^(*--' \\V''+'giu){r)\\dT I',^, + 



Here, for I^ ^, since 
by Lemma 1 5. 11 one has 



' (l + i)-i-4ln(l + t), n = 2, 



< C X 



+ ?^>3. 



For /g 2 , similar to (j4.5p , one has 

7^ 2 < C y^-^(*-^)||w(T)||Loo||V'=+iu(r)|ldr < + 

By putting estimates on I'^ ^ and /g 2 into (j4.7p , one has 
(4.8) 4<Cp(0(l + i)-^-'^. 

Thus, (|4.2p together with (|4.6p and (|4.8p yields (|4.ip and hence completes the proof 
of Lemma 14.11 □ 



Now, we turn to the case when s > 1. To the end, similar to define iVi in p.25p . 
let us define the integer N2 in terms of n > 2, s > 1 and N as follows. Denote 

71 + 1 

(4.9) li{n, s,k)^k + l + [(^— + fc + l)(s - 1)] + , 

r (fc + l)([s]+-l), n = 2, 

(4.10) 4(s,fc) =4(s,fc) + i/(s,fc), v{s,k):=\k{[s\+~\), n = 3,4, 

[0, n>5. 

Then, N2 is defined by 

0<fc< [[i^], 4(n,s,fc) <7V, ] 



(4.11) N2 = sup <^ fc e Z 



and £5(3, fc) < N 



Lemma 4.2 (case for s > 1). Lei 71 > 2, s > 1 a77.rf N > No, where Nq is given 
in (|3.24p . Suppose that ||Mo|!/fNnLi sufficiently small, andu{x,t) is a solution to 
the Cauchy problem (|l.ip obtained in Theorem ] Then, for any t >0, 

(4.12) ||v'=(u-e^^=*wo)(0|| <Cp(t)(l+t)-t-^, 

where < k < N2, and pit) is defined in the same way as in Lemma\4-.1\ 
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Proof. Take 1 < fc < iV2. Similar to obtain (|4.2p . one has 

(4.13) WV" {u{t) - e^P^'uo)\\ <h+h, 

where as for and /g in (|4.2p . 14 and also have the same definitions as I4 
and /5 given in p.35p . Once again, in what follows we would improve the previous 
estimates on I4 , /g with the help of the obtained time-decay estimate (|1.5p in the 
case when s > 1. For that is I4 in p.BSp . one can apply (|2.4p with k replaced 
by fc + 1 and with p = 1, r = g = 2 and = !i±iM ^ fe±i foj. n > 2, so that it 

follows 

(4.14) l';<C /'(l + i-r)-t-^||5(w)(r)||^,dr 

Jo 

+ C [\l + t- r)-^-'^ ||V^+i+[(^+'=+i)(^-i)l+g(ii)(r)||rfr 

:= I41 + /4 2- 

For the term I4 i, as in p.38p , one has 



(4.15) I'l, < C£°^P^{t) 



(1 + <)-!- 2 ln(l + i), n = 2, 
il + t)-f-^, n>3 

<Cf°P(Op(t)(l+i)-^' 



where the last inequality holds true since n > 2. For the term 142, by applying 
Lemma |5. 11 one has 



(4.16) l'l2<C' / '(l + t-T)-'^-'*^||u(r)|U=o||V^^("^^''=)«(T)||dr 



<CMo{t){l+t)-^-^ I {l + T)-^\\u{T)\\HNdT 



< CMo(t) V^jv(t)(l + t)-^^-^ I ' {1 + t)-^+Ut 

< CA/o(<)V^Ar(t)(l + t)- 



71 k + 1 
4 2 



where the definitions ([3^ . and (|49l) for A'/o(0: ^Nit) and 4(?i, s, fc) were 

used and we also used p.40p for all n > 2. Then, combining (|4.14p . (|4.15p with 
(|4.16p gives 



(4.17) /; < Cp(t)(l + t)-t-^ ^£°P^(t) + Mo{t)V£^ 

Next, we turn to estimate /g in (|4.13p or equivalently I5 in p.35p . As in p.42p . 
we write for simplicity that /g is bounded by the sum of /g ^ and /g 2 which have 
the same definitions as /s^i and 15,2, respectively. Then, similar to obtain p.43p . it 
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follows that 

f + ln(l + i), n = 2, 

(4.18) /5,i<CC^(t)x <^ 

[(l + t)-t-2, n>3 

where n > 2 was used for the last inequality. To estimate 2 ; notice that from the 
definitions (|4.10p and (|4.1ip of (^{Sjk) and N2, one can take 6 > small enough 
such that 

1 + [2(1 + S){s - 1)]+ + k[s]+ < 4(s, k)<N 
for any < fc < A^2- Then, it follows from Lemma [CTl that 

[|yi+fc+[2(i+5)(.-i)]+^(^)(^)|j < c\\u{t)\\l^ ||Vi+'=+P(i+^)(^-i)l+«(r)|| 

{\\V''+^u{T)\\^N-ik+i)M^, n = 2, 
||V^^i(T)||^„_M=i+, n^3A, 
l|V''~^M(T)||^iv-(fc-i)w+ , n > 5 

< CA/o(t)yMO(l + r)-5"'^+^ 

for any < t < t, where ri{k) = k + 1 for n = 2, k for n = 3, 4 and fc — 1 for n > 5. 
This further implies 



2 +4c(t- 



/5.2<C^^4Wv/^^ / (l+t-r)-(i+*)(l + T)- 



<CAfo(i)V^Jv(i)(l+^)~^ 



fc + l n + 2T|(fc)-2(fc + l)-l / 2 



■(l+T)-(l+*)dT 



<CMo(t)VMi)(l + 0~^' 



where n > 2 was used. Therefore, by combining estimates on /g i and /g 2 above, 
one has 

(4.19) i; < I'l, + i;., < Cp{t)(l + tr'i-"-^ [£°l (t) + Mo{t)^£^) . 

Thus, together with (|IT7)) and yield that for any < t < T, 

II V''^ {u{t) - e^P'*uo) \\<C [f]^P(t) + A/o(t) VM^)] P{t){l +tr^-'^ 

< Cp(t)(l+t)~T-^ 

with < k < N2, where the uniform-in-time boundedness of £'^^{t), MQ(t) and 
£N{t) was used. This proves (|4.12p and hence completes the proof of Lemma 

□ 



Next, for the initial data uq given above, we define the desired time-asymptotic 
profile u* ~ u*(x,t) by 

u* {x,t) = G{x,t + 1) / UQ{x)dx, 

where G = G{x,t) = (i^snt)^^ e is the usual Green function of the linear 
heat equation. As used in [52], we have the following well-known result. 
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Lemma 4.3. Let n > 1, k > and 1 < q < 2. //0 G L"? (R"), then 

||V^G*0(i)|| < Ct-*(i-^-*)||0|U, 
for any t > 0. Also, if 4> E L\ (M") and Jjg„ (l){x)dx = 0, t/ien 

||V^G*0(i)|| <Ct-5(i-^-'^)||0||i;, 

for any t > Q. 

Based on Lemma lOl one can show that e'^='^*uo is well approximated by u*{x, t) 
time-asymptotically. 

Lemma 4.4. Let n > 1, k > 0, and let uq G L\ (M") and Jjj„ uo{x)dx ~ 0. Then, 

(4.20) II V*^ (e^^^*uo - u*) (i)|| < C(l + i)-^-^ 
for any t > 0. 

Proof. Define 

(t)o{x) = G{x,l) = {4h,7t) 2e ^^-^ . 
Then, G{x,t + 1) = G{t) * (j>o{x) holds true. Consequently, one can write 

e'''^\io-u* = G{t)*{uo~ M(j>o), M := uo{x)dx. 

Note that /ju„ 4>Q{x)dx = 1 and /][{„(uo — M(f>Q)dx = 0. Therefore, for fc > 0, by 
applying Lemma [131 we deduce that 

IIV*^ (e'^^^^uo -u*(t))|| = ||v'=G*(uo-A/0o)(t)|| <G(l + 0"^"^ 
for any t > 0. This proves (|4.20p and hence completes the proof of Lemma l44l □ 

Proof of Theorem II. 3t For the solution u to the Cauchy problem and the 
desired time-asymptotic profile u*, their difference can be rewritten as 

(4.21) u- u* ^ {u- e^P^'uo) + (e^^=*Uo - e'^=^*Uo) + (e^-'^*Mo - u*) . 

Suppose that ||Mo||/f]VQ^i is small enough and uq G L\ with Jg„ u^dx ~ 0. When 
s < 1, combining of (piU)) . with (|i:iT|) yields ([Til) for < A: < iV - L 

When s > 1, combining of fiT^ . ^(TH) . with leads to (HH]) for 

< k < N2- The proof of Theorem 11.31 is complete. 

5. Appendix 

In this section, we prove some inequalities about L^ upper bounds of some non- 
linear terms, which have been used in the previous sections. This first inequality is 
about the L^ estimate on any two product terms with the sum of the order of their 
derivatives equal to a given integer. 

Lemma 5.1. Let n > 1. Let ~ (a}, ■ • ■ , a^^) and c? = (a^, • ■ • , a^) he two 

multi-indices with \a^\ — fci, |q;^| — k2 and set k = ki + k2. Let 1 < p, q,r < oo 
with 1/p =l/q + l/r. Then, for Uj : R" R (j = 1, 2), one has 



(5.1) 



LP(K") 



< C-l hlllL^CR") lr^2||^wR„, + h2|li,(H.) ||V"^Ul||^. 



for a constant G independent of ui and U2. 
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Proof. It is similar to the proof for the case of n = 1 given in (TOl Lemma 4.1]. 
Here, for the convenience of readers, we present a complete proof for the multi- 
dimensional version. First, for j = 1, 2, set 9j = ^ and define pj by 



1 _ 1 - 
Pj 



q 



r 



It is straightforward to verify ^ = ^ + ^ since 61+62 = 1. Therefore, applying 
the Holder inequality, the Gagliardo-Nirenberg inequality and the Young inequality, 
one has 



||5"^MlC'"^U2| 



< 








9"^2 










LPi(I 


t") 






5") 



= C (hl|L,(K.) II V'=U2||^,(j,„)) ' [\\u2\\mm-) ||VV|L.(K„) 



IV^Mll 



< (hl|lL<,(R") ||V U2||^,(^„^ + ||lt2||i,( 

where the following Gagliardo-Nirenberg inequality over R" was used: 



|V''^u,|Lp,,„„, < ||v''wj||/,,^„. \\u 



"■JUL"] (R") 



) II"JIIl<!(K") 



with ± ^ + 6, (i - 
Lemma [521 is complete. 



(1 - 6j)^, j = 1,2. Then ([SH) follows. The proof of 



□ 



The second inequality is the generalization of Lemma 15.11 up to the case of 
products of several terms. 

Lemma 5.2. Let n > 1, I > 2 be integers. Let — (aj, ■ • • I < j < I, be 

multi-indices with |q;-' | = fc-,-,1 < j < ^, and set k = ki + - ■ ■ + ki. Let 1 < p, g, r < 00 



with 1/p = 1/q + 1/r. Then, for u = (ui, • • • , ui) 



has 



(5.2) 



<c^I|u||Um")II"IIl'(M")II^'"I 



LP( 



for a constant C independent of u. 

Proof. Once again, it is similar to the proof for the case of n = 1 as in [131 Lemma 
4.1]. In fact, for j = 1, • • ■ , set 6j ~ ^ and define pj by 
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Notice ^ = ^ + ' ■ ' + ^ since 9i + ■ ■ ■ + di = 1. Therefore, applying the Holder 
inequality and the Gagliardo-Nirenberg inequality, one has 



< 



<cn(iiu,i 
<cn(ii"i 



II "JIIl-cr") 



L°°(R") II'-''JI1l9(R") II ^ "j|Il'-(R") 



L'-(R") 



|V*-u|| 



lL'-(R") ' 

where similarly before, we used the Gagliardo-Nirenberg inequality over . 



I V'^J?; 11 < llv'^7; " ^ 

I ^ "jIIl''3(r") ^ Ir " 



'J|Il'-(R") ii"jIIl<!(R") 



1 /s, 



1 fc 



with 



for each j = 1, • • • , / . This proves Lemma 15.21 



+ +(1-^^) 



1 



a -1)9 



□ 



From Lemma [5.2[ we have the following corollary, which provides a limit situation 
of (lOI). 



Corollary 5.1. Let n > 1 and 1 < p < oo. Let a = (ai, • • ■ , a„) 6e a multi-index 
with \a\ = A:. Assume that F{u) is a smooth function ofu. Then, there is a constant 
C(||u||l'^(K")) depending only on with C(e) when e — >■ smc/i f/iai 

\\d" (f(kKJ||^,(r„) < c(hiu^(R„)) ||v'=+i7.||^,(j,„) 

for all 1 < i < n. 

Proof. By Leibniz formula, for each 1 < « < n, 

1/^1=1 

fe l^ai I 



l/3| = l '=1 q1 + -+q'=/3 J = 1 

|ajl>l,l<i<' 

Define a'+i = a - + Ci = a - + (0, ■ • ■ , 0, 1, 0, ■ • • , 0). Then, it follows from 
Lemma 15.21 that 



w (^^(«k.)iIlp(r.) < \\F{u)^"u,^\ 



k \0\ 
|/3| = 1 Z=l 



i+1 



OPTIMAL DECAY RATES TO SOME CONSERVATION LAWS 
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This proves Corollary 15. II □ 
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